Temporal Distinguishability of an N-Photon State 
and Its Characterization by Quantum Interference 
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We present a multi-mode model to describe an arbitrary N-photon state with a wide spectral range 
and some arbitrary temporal distribution. In general, some of the N photons are spread out in time 
while other may overlap and become indistinguishable. From this model, we find that the temporal 
(in) distinguishability of photons is related to the exchange symmetry of the multi-photon wave 
function. We find that simple multi-photon detection scheme gives rise to a more general photon 
bunching effect with the famous two-photon effect as a special case. We then send this N-photon 
state into a recently discovered multi-photon interference scheme. We calculate the visibility of the 
multi-photon interference scheme and find that it is related to the temporal distinguishability of the 
TV photons. Maximum visibility of one is achieved for the indistinguishable N-photon state whereas 
the visibility degrades when some of the photons are separated and become distinguishable. Thus 
we can identify an experimentally measurable quantity that may quantitatively define the degree 
of indistinguishability of an N-photon state. This presents a quantitative demonstration of the 
complementary principle of quantum interference. 

PACS numbers: 42.50.Dv, 03.67.Mn, 42.50.St 



I. INTRODUCTION 

The coherence properties of an optical field are best 
described by the field correlation function in space and 
time Most commonly used quantity to characterize 
the coherence property of an optical held is the coher- 
ence time or coherence length for temporal coherence. 
Roughly speaking, the coherence length of an optical held 
is the distance within which the held can be described as 
a single uninterrupted wave train. In other words, any 
two points within the coherence length will have a hxed 
phase relationship. However, this description is primar- 
ily concerned with wave aspect of an optical held and 
is based on the interference effect observed in intensity 
or single photon interference effect. More specifically in 
terms of the quantum coherence theory it is related 
to the held correlation function of 



r( T ) = (£(-)( t + T )£(+)(i)) 1 



where 



£(+)(*) 



>2it 



duja(uj)e 



(1) 



(2) 



for a quasi-monochromatic held [3| and the average is 
over the quantum state of the held. The visibility of the 
single-photon interference fringes is simply the absolute 
value of the normalized held correlation function [i[ : 



7 (r) = r(r)/r(o). 



(3) 



However, this description becomes rudimentary when we 
start to deal with the cases involving more than one pho- 
ton in quantum information. One may use a higher or- 
der correlation function such as the intensity correlation 



function 0] 

= {M-\ti)...M-\t N )M+\t N )...E^(t 1 )). (4) 

which is related to an N-photon coincidence measure- 
ment. However, this function does not provide any in- 
formation about photon entanglement, i.e., quantum su- 
perposition of different states. 

The realization of multi-particle entanglement is 
paramount in achieving most of the tasks in quantum 
computing and quantum information processing 0, 0|. 
While there are many ways to create entangled multi- 
particle state, the straightforward method is to start from 
independent single photons 0. Knill, Lahamme, and 
Milburn Q have shown that quantum computing can be 
realized with single photons and some linear optical ele- 
ments via multi-photon interference. N-photon entangle- 
ment is thus produced from single-photon states. This is 
one of the primary reasons behind the big rush in creating 
light sources with single-photon on demand 0, [l(| ■ 

While most of the analysis are based on the single mode 
model, i.e., all the photons in one single temporal mode, 
this is, on the other hand, impossible to achieve in exper- 
iment. The multi-mode nature of light inevitably reduces 
the effect of photon interference and leads to degradation 
in information processing. One often uses the fidelity 
quantity of quantum states to characterize the degrada- 
tion. But this description has emphasis only on the end 
result of the process and spares the true culprit of the 
process, that is, the multi-mode nature of light. 

For monochromatic held of only one frequency com- 
ponent, the held can be represented by an infinite wave 
train. Photons can appear anywhere in this wave train 
and are indistinguishable from each other. They will pro- 
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duce maximum effect of entanglement. However, when 
many frequency components are excited, an optical field 
is no longer monochromatic and the wave train becomes 
finite with a length of the order of the coherence length 
of the field. With multiple photons, we generally cannot 
use a single wave packet to describe them. We cannot 
assign separate wave packets to describe each photon, ei- 
ther. This is because of the possibility of multi-photon 
entanglement. Thus, an issue is raised about how to de- 
scribe the different situations of temporal distribution of 
photons and distinguish these situations experimentally. 
Recently, this issue was addressed in the four-photon 



case [111 [12j for distinguishing a genuine four-photon 
polarization entangled state from a state made of two 
well separated pairs of photons. The difference lies in 
the multi-photon interference: an entangled four-photon 
state will give rise to the strongest multi-photon interfer- 
ence effect whereas two well separate pairs produce less 
interference effect. This is in consistence with the com- 
plementary principle of quantum mechanics which states 
that quantum interference is a result of indistinguishabil- 
ity of the paths but if the the paths are distinguishable, 
the interference effect will be gone. Partial distinguisha- 
bility will lead to reduced interference effect, as described 
by Eq.([3]) in the coherence theory for the single-photon 
interference. Four-photon interference experiments were 
performed to distinguish an entangled four-photon state 
from two independent pairs of photons [III, UH, 0, H3 | 

However, the above mentioned interference scheme on 
the four-photon state cannot be generalized to arbitrary 
photon number. More recently, Sun et al [HI, [l6| and 
Resch et al (l7j independently constructed a quantum 
state projection measurement scheme and applied it to 
a maximally entangled N-photon state (the so-called 
NOON state) for the demonstration of multi-photon dc 
Broglie wavelength without a NOON state, ft turns out 
that this new projection measurement scheme is based 
on a multi-photon interference effect that depends on the 
temporal distribution of the photons involved. Since the 
new scheme can be easily generalized to arbitrary photon 
number, it can be used to study the relation between the 
multi-photon interference effect and the temporal distin- 
guishability of an N-photon state. We will show that 
the various scenarios of temporal distribution of photons 
give rise to different visibility in the multi-photon inter- 
ference, which provides a direct measure of the degree 
of temporal distinguishability of a multi-photon state in 
a similar fashion to the coherence theory [Eq.Q]. This 
is a quantitative investigation into the complementary 
principle of quantum interference. 

In the following, we will first review the two-photon 
and four-photon cases to look for the relation between 
temporal distinguishability and multi-photon interfer- 
ence. We then will generalize to an arbitrary N-photon 
state and present the criteria for photon indistinguisha- 
bility and distinguishability. In Sect. IV, we use quan- 
tum coherence theory to calculate the result from a di- 
rect N-photon coincidence measurement and discuss the 



generalized photon bunching effect. This measurement 
process is not sensitive to the different temporal distri- 
bution of the photons. In Sect.V, we introduce the newly 
constructed NOON state projection measurement and 
demonstrate how it can be used to characterize the de- 
gree of temporal indistinguishability for the simple three- 
photon case. We will generalize the discussion for three- 
photon case to the general N+l-photon case. In Sect. VI, 
we will discuss an even more general case and present 
the numerical results for a few special cases. We con- 
clude with a discussion. 



II. TEMPORAL DISTINGUISHABILITY FOR 
THE CASE OF TWO PHOTONS AND FOR THE 
CASE OF TWO-PAIRS OF PHOTONS 

The first discussion about the t emp oral distinguisha- 
bility was by Grice and Walmsley [18| , who investigated 
the visibility in a Hong-Ou-Mandel interferometer [l9j] 
with a two-photon state input from type-II parametric 
down-conversion. Later on, Atature et al [2(j performed 
experiment and confirmed the degradation of the two- 
photon interference visibility predicted in Ref. [lq due to 
temporal distinguishability. 

In the discussion of Ref. [l8| . the multi-mode descrip- 
tion of the two-photon state is given by 



$2) = / <kjxdw2$(oJi,(J2)al(uJi)al(uJ2)\0) > 



(•5) 



where s, i denote the two correlated signal and idler pho- 
tons from parametric down-conversion. For type-II pro- 
cess, we have $(it;i,u; 2 ) 7^ ^{^2, ^>i) due the birefringent 
effect of the nonlinear crystal on the ordinary and extra- 
ordinary rays. The maximum visibility in the two-photon 
Hong-Ou-Mandel interferometer has the form of 

V9 = M = / ^1^2$* (^1, ^2)^(^2, uji) 

Jrfwidw 2 |*(wi,w 2 )| 2 ' ' " 

M 2 is defined as a degree of permutation symmetry. Note 
that M 2 = M| and < \M 2 \ < 1. The visibility or the 
degree of permutation symmetry is one if and only if 
LU2) satisfies the permutation symmetry relation: 



(7) 



As stated in Ref. [181 ] , this permutation relation is a signa- 
ture of spectral indistinguishability of the two photons, 
that is, we cannot tell the difference between the two pho- 
tons through their spectra. This in turn gives temporal 
indistinguishability if we consider the Fourier transfor- 
mation: 



G(*i,i 2 ) = J daj 1 duj 2 ${uj 1 ,uj 2 )e~ 



1(^1*1+1^2*2) 



(8) 



Combination of Eqs.JTl [8]) gives directly the symmetric 
relation: 



G(ti,t 2 ) = G(*2,ti), 



(9) 
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for all times of t\, t 2 . 

On the other hand, the visibility is zero if $(w 1; W2) 
does not have any overlap with $(w 2 ,wi), which is char- 
acterized by the orthogonal relation: 



du>idu>2®* (^1, U) 2 )$>(ui2, tvx) = (10) 



or in time 

J dt 1 dt 2 G*(t 1 ,t 2 )G{t 2 ,t 1 ) = 0. (11) 

This orthogonal relation indicates that the two functions 
G(t\,t 2 ), G{t 2 ,t\) have no overlap. 

At this point, it is not easy to see what is the physical 
meaning of Eq.JTl]). However, if we go back to Eq.([7]) 
and introduce a non-symmetric factor of e luJ2T , we find 
that the equivalent $ (0^1,0^2) in Eq.© in this case will 
be (l>i,u) 2 ) = $(oji, L02)e lU2T : which is not symmet- 
ric with respect to u)x,u> 2 even if $(wi,tt>2) is symmetric. 
This extra phase can be introduced by acting the evo- 
lution operator U(T) — exp(— iu)2d\diT) on the state in 
Eq.© for an extra free propagation time T of the idler 
photon. This then creates a time delay T between the 
two photons. Then the visibility in Eq.© becomes 

vat) - S^^^u^m^^y^-^ 

J duJiduj2\<&(uJi,uJ2)\ 2 

Notice that if the delay is large enough [T >> T c ~ 
1/Aujpoc with Awpcc as the range of <f>(u)-L,u> 2 )], we 
will have V 2 (T) = or &(wx,U2) satisfies Eq.([T0"|). Since 
T is the relative delay between the two photons be- 
fore they meet at the beam splitter of the Hong-Ou- 
Mandel interferometer, we may believe that there is a 
large enough delay between the two photons so that the 
two photons become distinguishable in time when they 
arrive at the beam splitter. So the orthogonal relation 
in Eq. (fTT)| or Eq. fTTj) corresponds to the situation when 
the two photons are well separated in time and form two 
non-overlapping and distinguishable wave packets. 



Therefore, the visibility in the Hong-Ou-Mandel in- 
terferometer in Eq.® is a direct measure of temporal 
distinguishability of the two photons. This is very much 
similar to the role of the field correlation function 7 of 
Eq.([3]) in defining optical coherence of a field. 

For the four-photon case, temporal distinguishability 
between two pairs of photons was first studied by Ou, 
Rhee and Wang [HQ in a similar scheme as the Hong- 
Ou-Mandel interferometer but with four photons. It was 
found that the visibility in four-photon interference is 
directly related to a quantity £ /A, which is a measure 
of the temporal distinguishability of photon pairs from 
parametric down-conversion: when £ /A << 1, the pairs 
are well separated from each other corresponding to the 
so-called 2x2 case but when £/A = 1, the two pairs 
are overlap in time and form an indistinguishable four- 
photon state corresponding to the 4x1 case. 

From the definition of the quantities £ and A in 
Ref . [l4| . we rewrite them as 

£ = y dujid^duj'iduj^* {tui, u>2)^*{^[, u' 2 ) 

x$(wJ,W2)S(wi,w£), (13) 

and 

A = J duj 1 duJ2did[did 2 \$(id 1 ,uj2)$(uj' l7 uj' 2 )\ 2 , (14) 

where <f>(u>x,cj 2 ) is the two-photon wave function in 
Eq.©. 

On the other hand, the four-photon state from Ref. [3] 
has the form of 

l^) = J duJiduJ2duj' 1 dLd 2 $>4(uJi,UJ2',Lij' 1 ,<JJ 2 ) 

xat( Wl )al(w2)atK)«!K)|0>, (15) 

where $4(0;!, u> 2 ; u>' 2 ) = w 2 )$(cc4 , w 2 ). Then we 

can rewrite the expression for £ and A in Eqs.(fT3l [M]) 
and obtain the quantity £ /A as 



£ J d(jJid(jJ2dLj[dlJ 2 $l(ui,UJ2,u[,UJ 2 )&4(lj[,U2;Ui,U 2 ) 

A J duJi<faj 2 <faj' 1 du> 2 \<S>4(u)i, cl> 2 ; oj[, (^> 2 )\ 2 

I 



(16) 



Recall that this quantity is a measure of the temporal dis- 
tinguishability of two pairs of photons. But from Eq. (|16[) . 
we find that this quantity is again dependent on the per- 
mutation of the wave function $4(^1, 0J2] lu[, w' 2 ) similar 
to that in Eq. © and it is one if and only if we have the 
permutation symmetry of 



Therefore from the discussion on the meaning of the 
quantity £ / A, we find that the symmetry relation in 
Eq. (fT7j) corresponds to the case when the two pairs are 
completely overlap in time and become temporally in- 
distinguishable (the 4x1 case) whereas the orthogonal 
relation 



$4(^1, W2, u[, U>' 2 ) = U 2 \ Wl, U! 2 ). 



(17) 



dLUldLU2duj' 1 duJ 2 ^>l(LUl, u> 2 , U^UJ^) 
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x$ 4 K,w 2 ;wi,w 2 ) = (18) 

leads to the case of completely separated pairs of photons 
(the 2x2 case). 

From the experiments and analysis on four-photon in- 
terference with two pairs of photons by parametric down- 
conversion [TH, [TJ, [l5| , we find that the visibility is not 
zero even for £ /A = 0. This can be attributed to the ex- 
istence of two-photon interference since we usually have 
two-photon indistinguishability with exchange symmetry 
in Eq. j2J . Note that £ j A concerns the permutation sym- 
metry between two different pairs, i.e., exchange between 
the group of {wi,^} and the group of {u>[, uj' 2 }- The ex- 
change within each group is symmetric due to Eq. ([7]) . 

Next we will generalize Eqs.(0 QU) and Eqs.jTUl QH) 
to an arbitrary iV-photon case and relate them to the 
visibility of some iV-photon interference experiment. 

III. DESCRIPTION OF A TEMPORALLY 
DISTRIBUTED N-PHOTON STATE 

Now we can generalize Eqs.Q HUf of the two-photon 
case and Eqs. (fT7l IT8]) of the two-pair case to arbitrary N 
case. An arbitrary N-photon state of wide spectral range 
can be generally described by 

\$n) = Af~ 1/2 J dujiduj 2 ...dujN${ui, wjv)x 

xa t (tu 1 )a t (uj 2 )...a ji (Lu N )\0), (19) 
where the normalization factor Af is given by 

N = J dldidld2---duJN&* (idi, 07jv)x 

x^$(P{ Wl ,.., WN }), (20) 
p 

where P is the permutation operator on the indices of 
1,2,..., N. and the sum is over all possible permutation. 
There are totally N\ terms. So the value of Af ranges 
from I to N\I with I = j du>xdu>2- . .du> n\$(uji, ...,u}n)\ 2 . 
The maximum value of N II is reached when 

$(wi,...,Wit) = $(P{wi,..,wjv}) (21) 

for all P. Similar to Eqs. (0 [T7|) . this corresponds to a 
case when the N photons are indistinguishable in time. 
We refer to this case as the JVx 1 case, meaning that all 
N photons are in one indistinguishable temporal mode. 
This single-mode description of an N-photon state is 
more vivid in the special case when ^(wi, ...,wjv) is fac- 
torized as 0(wi)(/>(w2)...0(wjv) and the N-photon state 
simply becomes 

l^) = ^(0) tAr |O)-|iV) (22) 

with 

A(<(>) = I du(t>{u)a{w) ( j dio\(t)(uj)\ 2 = 1). (23) 



Note that A(<p) satisfies [A, A 1 ] = 1 and represents the 
annihilation operator for a single temporal mode char- 
acterized by <j){ijj). The single-photon state has a 
single-photon detection probability of |(?(t)| 2 with a tem- 
poral shape of 

g(r) = -L / d^{oj)e-^ (24) 

V27T J 

and normalization relation 

J dr\g(r)\ 2 = 1. (25) 

The other extreme case of Af = / requires 
$(wi, u>jy) be orthogonal to all the permuted functions 
$(P{wi, ...,Wjv}) in the similar ways in Eqs. fTUl [Tg|) and 
thus corresponds to the situation when all photons are 
well separated in time. We refer to this case as the 1 x N 
case, meaning that each photon is in its separate tempo- 
ral mode and there are totally N independent modes. 

For the situations in between the two extreme cases, 
the value of Af is between I and Nil. For example, as- 
sume that the spectral amplitude $({w}) have partial 
permutation symmetry, that is, 

$(wi,...,wat) = $(P {ni} {wx,...,u N }), (26) 

where the permutation P{ ni } only applies to a subgroup 
of {lo\,uj2, ■■■,^n}- In the meantime, it also satisfies the 
orthogonal relations: 

J dw 1 ...du N & : (ui, ...,uj N )^>(P lj {uj 1 , ...,uj n }) = (27) 

for permutation P^ between different subgroups ({rii} 
and {nj} 7 i ^ j) defined in Eq. ([2"6"]) . Then it can be 
easily shown that Af = ni!ri2!...nfc!P In the simple case 
when $(w) can be factorized as 

$(W1, U N ) = 0i(wi)...0i(cJ„ 1 )(/)2(w„ 1+ i)x 

X...^ 2 (^n 1 +n 2 )—4>k(^N) (28) 

with the orthogonal relations 

f du^tffaW^M^M 1 ^) = ( l ^ i). ( 29 ) 

the N-photon state in Eg. ([IT?]) becomes 

\®n) = —.\ni)^—An 2 )fc...^—\n k )<i, k . (30) 
n\\ n 2 \ nkl 

This is the situation when the N photons are divided into 
k subgroups with m(i — 1, 2, fc) photons in each group 
in a single temporal mode characterized by <pi. This sit- 
uation is denoted as n\ + ... + case. 

For simplicity of later argument, let us consider an- 
other special kind of N-photon state with 

= ^i)e ,UlTl ...^K)e w " T ". (31) 
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With this <&, the N-photon state can be viewed as direct 
product of N identical single photon wave packets: 



\N) T = |Ti)® |T 2 )®...® ITjv), 



with 



\Ti) = I dw0(w)e iwT 'a t (w)|O) 



(32) 



(33) 



This state can be viewed as from single-photon sources 
such as quantum dots (see below for details). However, 
the quantum state in Eq. (|3"2"|) is not normalized. Substi- 
tuting Eq. (|3Tj) into Eq. (|20|) . we have the normalization 
factor as 



fe 



(34) 



When T% = T% = ... = TV, we recover the case when all 
N photons are in one single temporal mode with Af = N\ 
(Nxl). On the other hand, if \T t -T j \» 1/Au(i ^ j) 
with Alo as the bandwidth of 4>{lu), we have Af = 1. This 
is the case when all the photons are well separated from 
each other (1 x N). 

The N-photon state in Eq. (fT5)) describes a state when 
all photons are in one spatial and polarization mode. 
They only differ in spectral mode. In practice, although 
N = 2 case can be easily obtained from degenerate para- 
metric down-conversion, such a state with N > 2 is not 
easy to produce directly. It can be produced indirectly 
from single-photon states with a set of beam splitters as 
shown in Fig.l, where the single-photon sources are, for 
example, quantum dots. The quantum state for the input 
fields has the general form of 

= j ' duj 1 duj2---duJ N <^{uj 1 , ...,cu N )x 

xd\(cj 1 )dl(uj2)---^ N (LJ N )\0) (35) 



with the normalization relation: 



dtuiduj2---dLL)pj(fr* (u>x, ...,u>n)&{ui, ■~,Wn) = 1- (36) 



Here d^(j = 1,...,N) is the creation operator for each 
input mode. Photons are possible to exit at any of the N 
output ports. To produce a state of the form in Eq. lfTT))) . 
however, we only consider the possibility when all N pho- 
ton exit at one port, say, 6/v port. It is straightforward 
using the beam splitter theory to show that the projected 
state is 



I* 



1 



N N / 2 



dujiduj2---dujN^(y>i, ...,Wjv')x 



x&yw 1 )&yw a )...&k(w JV )|o>, (37) 



N out — 



which is in the form of Eq. (|T9|) . This state is not normal- 
ized because it is a projected state with the probability 
of projection as PQ$n}) = N )out\\ 2 = N /N N . The 
delay factors {e^} in Eq. (j3T]) can be easily introduced 
on individual mode dj before the beam splitters via the 
free-field evolution operator Uj(Tj) — exp(— iuijd^ajTj). 



/ / 
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a N 



i*„> 



FIG. 1: Generation of an N-photon state from single-photon 
states by beam splitters. 

More generally, to include different spatial and po- 
larization modes, the N-photon state has the following 
shape 



\<S> N )=JV- 1/2 J ^^..^V---^ (38) 
where {u^ 1 - 1 } = w x , . . . , Lo^h} , etc. The normalization factor Afk takes the form of 




Afk now ranges from / to n\\...nk\I. The special case when ^({o/ 1 )}, |q/ fc )}) factorizes is similar as before. 

IV. DIRECT N-PHOTON MEASUREMENT: herence theory in Eq.®. The average is over the quan- 

PHOTON BUNCHING EFFECT FOR N turn state of the system given in Eq.(Q2]| for an arbitrary 

PHOTONS N-photon state. For simplicity, we first apply it to the 



Next, we consider an N-photon joint measurement with 
the joint probability density given from the quantum co- 
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state in Eq. (j32j) . 

To carry out the quantum average, it is easier to first 
find the N-photon detection probability amplitude: 

CW(t 1 ,i 2 ,... ! t A r) = <0|i (+ )(t JV )...^ (+) (ti)|$iv).(40) 

Then rW(ti,t 2> ...,fw) = \C^(t u t 2 ,...,t N )\ 2 . From 
Eq.© for the field operator and Eq. ([5l]) for $, it is 
straightforward to obtain 

C w (ti,t 2 ,.,t w ) 

= Y,P[9(h-Ti)-9(t N -T N )], (41) 

p 

where the permutation operation P is on tit2-..tjsr and 
there are ./V! terms in the sum. 

The overall probability of detecting N photons 
together (N-photon coincidence) is an integral of 
rW(ii, t 2 , £jv) over all times ti, tjq: 

P N = I dt 1 ...dt N T {N \t u t 2l ... 1 t N ) 



= / dt x ...dt 



N 



J2P[9(h-T 1 )...g(t N -T N ) 



(42) 



J 



In the extreme case when T\ = T 2 = ... = TV, we 
obtain Pn(N x 1) = (AM) 2 T while in the other extreme 
case when \T t - Tj\ » 1/AO, we have P w (l x N) = 
AMT. Therefore, we seem to have 



or 



P N {N x 1) = AMPv(l x TV), 



Pv(AT x 1)/Pjv(1 xN) = N\, 



(43) 



(44) 



that is, the N-photon detection probability is N\ larger 
in the case of N identical photons than in the case of N 
separated photons. This can be thought of as the Bosonic 
photon bunching effect for N photons. The case of N = 2 
gives the familiar photon bunching factor of 2. 

However, as we know, the N-photon state in Eq. (|3"2"|) 
is not normalized. With the normalization factor consid- 
ered, we have instead 



P N (N x 1) = Pv(l x N) = AM. 



(45) 



For the case in between the two extreme cases, we may evaluate Eq. 



as 



Pv = f dt 1 ...dt N Y, p '[9*(h-T 1 )...g*(t N ~T N )]Y J P[9{h-T 1 )...g(t N -T N )]. 



(46) 



Since the sum is over all permutations, the integral does not change if we make the variable change: {ii...^} 
P'{t x ...t N }, i.e., 



p N = I dti...dt N g*{tx - T x )...g*{t N - T N ) ^ P[g{tx - r x )... fl (t w - T N )\ 
P , J P 

= AM f dh.-dtNg*^ ~ Tx)...g*(t N - T N ) £ P[g(h - T x )...g{t N - T N )}. 

J D 



It can be further shown that 



f dt 1 ...dt N g*(t 1 -T 1 )...g*(t N -T N )J2P[g{t 1 -T 1 )...g(t N -T N )} = J\f, 



(47) 



(48) 



where Af is given in Eq. (|34p . Thus we have 

P w = NIAf. (49) 

For a normalized N-photon state, we have Pjv = AM in 
any case as in Eq. ([45|) . 

For the multi-spatial and polarization state in Eq. (|38l) , 
we may find Pjv after some lengthy manipulation as that 
leads to Eq.gl]): 

P N = n 1 L..n k W k (50) 



for the un- normalized state and P 4 = nil.-.n^! for the 
normalized state. 

Hence, it is impossible to characterize different cases 
of temporal entanglement with just simple direct multi- 
photon detection for the normalized state. Furthermore, 
even for the un-normalized state, we cannot explore the 
temporal indistinguishability among different spatial and 
polarization modes with multi-photon detection, for J\fk 
depends only on the permutation symmetry within pho- 
tons in one spatial and polarization mode. 
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Before we proceed further, it is interesting to evaluate 
the multi-photon detection rates in some special cases. 
For example, for the single-photon detection rate Pi, we 
have 



P x = J dt($ N \&(t)E(t)\$ N ) 



(51) 



With some manipulation, it can be shown that Pi = NJV 
for an un-normalized N-photon state and P\ = N for a 
normalized N-photon state. 

The reason that we still discuss the un-normalized case 
of an N-photon state is because we encounter this kind 
of state in practice when a projection measurement is 
involved such as that in Fig.l. Consider, for example, 
a multi-photon state from degenerate parametric down- 
conversion, which, for small rj, has the form of [15] 



1$ 



PDC 



2 

) =C(\Q)+ V \<$> 2D ) + \\<$> iD ) +...), (52) 



with 



|<&2£>) = J duidu)2${wi, w 2 )a t (cJi)a t (cj 2 )|0) (53) 



and 



xa t (wi)a t (w 2 )a t (w' 1 )a t (w 2 )|0). (54) 

Here C in Ea.(|52p is a normalization factor but because 
\rj\ « 1, \C\ « 1 no matter what function $(oji,o; 2 ) 
is. Two-photon and four-photon detections project the 
state to ?7|$ 2 _d) and ?7 2 1 < &4Z?) /2, respectively, which are 
not normalized. From Eqs. (|2QI49[) , we then have 



P 2 = 2H 2 / dwidw 2 [|*(wi,w 2 )| ! 



+$*(wi,W2)*(w2,wi)]- (55) 

The last term is related to the permutation symmetry 
or the degree of two-photon temporal distinguishability 
and can be viewed as a two-photon bunching effect. For 
a state from parametric down-conversion in the degen- 
erate case as in Eq. (|52p . we usually have the symmetry 
$>(w\,u)2) — < &(o; 2 ,cj 1 ) so that 



P. 



2D 



4|?7| 2 / du)idu) 2 \${ui,U2)f 



(56) 



Similarly for four-photon case, we have 

P 4D = A8\i]\ 2 (A + 28) = 3P 2 (1 + 28/ A), (57) 

where 8, A are given in Eqs. (fT31 [T4")) . respectively. The 
dependence on 8/ A indicates that the extra term in 
Eq. (f57|) is a pair bunching effect - a generalized photon 
bunching effect for a multi-photon state. Direct mea- 
surement by Sun et al 16] confirmed the four-photon 
bunching effect in Eq. (|57|) . 



Another example is from non-degenerate parametric 
down-conversion in type-II x medium. The quantum 
state is similar to that in Eq. ([52|) 14{: 



2 

\1>NPDC) = C(|0) + V\$ > 2N) + Y I $ 4JV) 



(58) 



with 



and 



<S>2n) = j duj 1 duj2<f>(uj 1 ,Lj2)a 1 H (uj 1 )a\ / (LU2)\0) (59) 

| <&4jv) = J du)idu)2du)' 1 du)' 2 ^(u}i, w 2 )$(w' 1 , w 2 )x 

xaU^i)al(Lo 2 )a^ H (Lo[)aU^ 2 )\0)- (GO) 
From Eq. ([ST)]) . it is straightforward to have 

P2N = \v\ 2 J du>idLJ 2 Mu 1 ,u) 2 )\ 2 (61) 



and 



P 4N = 2\ v \ 2 (A + 8) =2P 2 2 (l + 8/A). (62) 



Although there is no photon bunching at two-photon de- 
tection, we still have the pair bunching effect that de- 
pends on the 8 j A quantity. 



V. N-PHOTON INTERFERENCE FROM AN 
N-PHOTON STATE 

As seen in the previous section, a direct N-photon de- 
tection scheme cannot characterize the temporal indis- 
tinguishability in the general case. Therefore, we need to 
seek another method. Since the direct result of photon 
indistinguishability is the interference effect, our scheme 
will be an N-photon interference scheme. As a matter 
of fact, a Hong-Ou-Mandel interferometer [l!| has al- 
ready been used to measure two-photon indistinguisha- 
bility from a type-II non-degenerate parametric down- 
conversion [l8|, |20| . Our method proposed in the fol- 
lowing will be a generalization of the Hong-Ou-Mandel 
interferometer from a two-photon case to an arbitrary 
N-photon case. 



A. NOON State Projection as a Measure for 
Distinguishability 

The NOON state projection measurement was recently 
proposed to demonstrate an N-photon de Broglie wave- 
length without the need for a NOON state It was 
demonstrated for N = 4 with states from parametric 
down-conversion [l6| and N — 6 for a coherent state [ItJ 
experimentally . The scheme is depicted in Fig. 2 where 
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the input is an arbitrary N-photon state of two polariza- 
tion modes in the form of 



N 

E 

fc=0 



c k \N-k,k). 



(63) 



The N-photon coincidence probability from the N detec- 
tors is proportional to 



P N oc \ (NOON\^ N ) 



(64) 



If the input state is of the form of 1 7V — fc, k)(k ^ 0, N), the 
output of the projection is zero. From the construction 
of the NOON state, we find this orthogonal projection is 
a result of N-photon interference and thus it can be used 
to characterize the temporal indistinguishability by the 
visibility in the interference. We will demonstrate this in 
the following sections. 



|9V> 



l 

AM 



FIG. 2: A NOON-state projection measurement. 8k = 2kn/N 
is the phase difference between H and V. bk oc Eh — Eve lSh . 



B. Three-photon case 

Let us start with a three-photon state of the form 
|2ff,ly). So the three-photon NOON state projection 
measurement should yield null three-photon coincidence 
in the ideal case when all three photons are in one tem- 
poral mode. However, there may be some delay between 
the vertical photon and the two horizontal photons due 
to birefringence. Furthermore, the two horizontal pho- 
tons may also be separated from each other. To account 
for the three scenarios described above, we cannot use 
the single-mode state of |2#, ly) and have to resort to 
the multi-mode model discussed in Sect. III. 

A multi-mode three-photon polarization state for 
|2fr, ly) has the form of 



|$3) = J d,LJldLJ2duJ 3 ^(iL!l,UJ2,LU 3 ) 



xa^(wi)a^(w 2 )a| / (w 3 )|0). (65) 

For simplicity of argument, we take u>2, w 3 ) in the 

form of Eq.jHT]): 



= 0(wi)e^ lTl 0(w 2 )e iW2T2 0(w 3 )e^ 3T3 . (66) 
We will use the un-normalized state because in practice, 
the state in Eq. ((65|) can be generated by superposing a 
weak coherent state \a) with a two-photon state \rj) = 
|0) + from non-degenerate parametric down- 

conversion: 



|*3> = \at)n\v) ~ |0)+a|l H ,0 v )+ 

+ (a 2 /V2)\2 H ,0 v ) + v \l H ,lv) + 
+ (a 3 /V6)\3 H ,0 v )+? 1 a\2 H ,l v ) 



(67) 



where the states are in a single temporal mode and we 
only write out states up to three photons. A three-photon 
coincidence measure like that in the N-photon NOON 
state projection will only have contributions from the last 
two terms. By making the coherent state weak enough 
so that \rj\ » \a\ 2 , we are left with only |2#, ly) term. 
Since \a\, << 1, the three-photon state is not normal- 
ized. 

For the scenarios presented in the beginning, we can 
relate them to different values of T u T 2 ,T 3 . So T X =T 2 = 
X3 is for the case of three photons all in one single tempo- 
ral mode. When |T 3 -Ti| >> 1/Auj,\T 3 -T 2 \ » 1/Aw, 
the V-photon is far from the two H-photons. When 
I jTJj — Ti\ >> 1/Acj, the two H-photons are far apart. 

For the projection measurement in Fig. 2 with N = 3, 
we have the electric field operators at three detectors as 



E (t) = [E H (t) 
Ei{t) = [E H (t) 
E 2 (t) = [Eh (t) 



E v (t)]/V6, 

e^^E v (t)]/V6, 

e^/ 3 E v (t)]/V6. 



(68) 



To find the three-photon coincidence probability, we first 
calculate the time correlation function 

r (3) (ii,*2,i 3 ) 

= {E^its)^)^!)^)^)^^))- (69) 

It is easy to calculate £^2(^i)£'i(^2)£ ; o(^3)| ( I > 3): 



^2(tl)£l(t2)£t)(t3)|*3 



—. ( E H E v E H e^/ 3 



E V E H E H + E H E H E v e^^ |$ 3 ). 



(70) 



6V6 V 

Here we dropped the terms that have no contribution. The order of the operators is kept for the time variables t^ti 
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With the state in Eq. ([65)) and $3 in Eq. (|66|) , it is straightforward to find 



^2(tl)4(*2)^(*3)|*3) = ^{ 



Gih^ta) + G(t 2 ,t u t 3 ) 



G(ti,*3,*a) + G(t 3 ,*i,«2) 



i2ir/3_ 



G(t 2 ,t3,tl)+G(t3,t 2) tl) H°> 



(71) 



where 



G(t u t 2 ,t 3 ) = — | ^ J du)idu) 2 du) 3 $(uJi,0J2,u 3 )x 



xe 



— i(u> 1 t 1 +U> 2 t2+U>3t3) 



= g(ti - T^gih - T 2 )g{t 3 - T 3 ) (72) 

with g(r) given in Eg . (f24[) . The three-photon joint detec- 
tion probability is an integral of the correlation function 
in Eq. ([69|) over all time variables: 



ft 



(73) 



We are now ready to discuss the three scenarios pre- 
sented in the beginning of this section. The interference 
effect is best measured by the visibility which is usually 
defined as the relative depth of modulation as compared 
to the situation when the interference effect is zero. In 
the three scenarios, we find the situation when the V- 
photon is far apart from the two H-photons corresponds 
to the case of no interference, which sets the reference 
line for evaluating the visibility defined by 



V, = 



\P 3 -P 3 (T 3 =00) I 
P 3 (T 3 =oo) 



(74) 



Experimentally, we can scan T3 from T 3 = 00 until we 
observe the dip in P 3 and use Eq. (|T4"|) to calculate the 
visibility. 

Depending on the separation between the two H- 
photons, we actually only have two distinct cases: (i) 
the two H-photons are completely indistinguishable with 
T\ = T 2 = T; (ii) the two H-photons are well separated 
and distinguishable in time with \T\ — T 2 \ >> 1 /Aw. 

In case (i) with T\ — T 2 = T, we have the exchange 
symmetry G(t 1 ,t 2 ,t 3 ) = G(t 2 ,tx,t 3 ) and Eq.$75§ be- 
comes after the time integral 



ft 



A 3 — £3 (AT) 



36 



(75) 



with AT — T 3 — T and 



A 3 = I dwidw 2 dw 3 \^(uJi 1 uj 2 ,ijj 3 )Y 
, 3 

|2 ' 



(76) 



£ 3 (r) 



duj\(j){<jj)\' 



dw\<j){ijj)\' 



(77) 



Note that £ 3 (0) = ^3 and £3(00) = 0. So from Eq.(j74]), 
we have the visibility for case (i) as 



V 3 (») 



|P 3 (AT = 0)-P3(AT = oo) 
ft(AT-oo) 



(78) 



The 100% visibility corresponds to the single-mode dis- 
cussion before. 

In case (ii) with \Ti — T 2 \ >> 1/Aw, there is no 
overlap between G(ti,t 2 ,t 3 ) and G(t 2 ,ti,t 3 ) so that 
/ dt 1 dt 2 dt 3 G*{t 1 ,t 2 ,t 3 )G{t 2 ,t 1 ,t 3 ) = 0. We obtain af- 
ter the time integral 



ft = 



^ 3 -£3(AT 1 )/2-£ 3 (Ar 2 )/2 
36 



(79) 



with Aft = T 3 — Tl and AT 2 = T3 — T 2 . So we will have 
two dips with half depth when T3 scans through ft and 
T 2 . The visibility of each dip is then 50%. 

In summary, we find that the scenario when the two 
H-photons are separated have a visibility of 50% while 
when the two H-photons are in one temporal mode, the 
interference visibility becomes 100%. Therefore, we can 
distinguish the two different scenarios in the three-photon 
case by measuring the visibility in the NOON state pro- 
jection measurement. Recent experiment by Liu et al plj 
realized the two scenarios described above and confirmed 
the corresponding visibility. Next, we will generalize this 
result to an N-photon state. 



C. N + 1-photon case 

Let us now generalize the conclusion in the previous 
section to the case of |1#, iVy) with an arbitrary integer 
N. The most general scenario in this case is when the 
single horizontal photon (H) is indistinguishable from m 
vertical photons (V) while other N — m V-photons are 
well separated in time from the m + 1 photons (the case 
of lHmV + (N- m)V or lHmV for short). The multi- 
mode description of this state has the form of 
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\$(lHmV)) = / (L>idL>2...du) N+ i$(u)i, ...,uj n ; uj N+1 )a t v (uj 1 )...a\ / (u) N )a ji H (uj N+1 )\vac) , 



with 



Here we take $ in the form of Eq. (|3Tj) for ease of calculation 



(80) 



(81) 



I 

When m H-photons are in the same temporal mode It is easy to first evaluate J?i(ii)...Pjv+i(ijv+i) 
with the V-photon, we have T\ = ... = T m = Tjv+i = T. \$(lHmV)) . After expanding the product, we find only 
But the other N — m V-photons are well separated from N + 1 nonzero terms of the form 
these m + 1 photons. This leads to \Tj - T k \ » 1/Au 
with j = 1, 2, ...,m, JV + 1 and k = m + 1, JV and the 
orthogonal relation: 



iV+1 



e iSk E v {t 1 )...E H {t k )...E v {t N+1 )\${lHmV)). (86) 



fe=i 

dhdhg^h-T^fc-n) 

Xg(ti - T k )g(t 2 - Tj) = 0. (82) For the state |$(LffmV)) in Eq.((80)). we have 



Now we are ready to evaluate the joint N + 1-photon 
probability Pn+i in the NOON-state projection mea- 
surement scheme with an input state of \<&(lHmV)) in 
Eq. (|80| . P/v+i is a time integral of the correlation func- 
tion from (N + 1) detectors: 

r w (ii,t 2 ,.,(«) 

= miHmV)\P N+i (t N+1 )...El(t 1 ) 

xE 1 (t 1 )...E N+1 (t N+1 )miHmV)), (83) 



E v (t 1 )...E H (t k )...E v (t N+1 )miHmV)} 

= G(Pk,N+i{h, ...,t N +i})\vac), (87) 



with 



S(ti, ...,*ivj ijv+i) — ^ G(P{<i, ijy}; t<v+i) 



with 



where 



and 



Ej(t) ocE v (t) -E H (t)e iS * + ..., 



E H ,v{t) 



du)aH,v{u)e 



(84) 



(85) 



J 



JV+1 



G(ti, ...,t N ;t N+1 ) = JJ g(t s -T s ), 



(89) 



where P^jv+i exchanges with ijV+i and P is a permutation of ti, tjy. For the case of WmH, we have 

m jV 



G(t 1 ,... ! ^ +1 ) = 5 (^ +1 -T)J] 5 (t s -T) J] ff(ti-T,) 



(90) 



s=l 



so that G(t\, t/v+i) bas exchange symmetry in t%, i m , £jv+i- The overall (JV + l)-photon coincidence probability 
is then given by 



P JV+1 (lJ?my) cx dt 1 ...dt N+1 e iS "G{P ktN+l {t u ...,t N+ i}) 
J k=l 

= Y^e* s *-*A J dt 1 ...dt N+1 g(p KN+1 {t 1 ,...,t N+1 })g*(p^ N+1 {t 1 ,...,t N+1 }). 

Diagonal terms of k = j in the double sum are all same because the integration is over all time variables: 
J dt 1 ...dt N+1 \Q(P kyN+ i{t 1 , ...,t N+1 })\ = J dti...dt N+ i\Q(ti, tjv+i) I • 



(91) 



(92) 
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Furthermore, 



dtx—dtN+i\G(ti, ...,tjv+i) 



dt\ ...dtiydtpj+i 



z2 p [9{h - Ti)...g(t N - T N )]g(t N+1 - T N+1 ] 



= f dt x ...dt N J2 P l9(ti ~ Ti)...g(t N - T N )} 



(93) 



where we used the normalization relation in Eq. p5| . From Eqs. (l4l)H4T)l) . we find that it is simply N\Af with AT given 
in Eq. (|34|) . So the diagonal terms of k — j in Eq. (|9Tj) are summed to be (N + l)N\Af. 
The cross terms in the double sum in Eq. ([9Tj) are given by 



e l{Sk ~ S]) I dt l ...dt N+1 g{p KN+1 {t u t N+1 \)g*{p^ N+1 {t u t N+1 }). 



(94) 



Let us consider one arbitrary term in the sum. The time integral part can be rewritten as 

^ / dti...dtN+iG(P{ti, tk-i, tN+i, tk+i, ijv}j tk) J3 G*(P'{ti, ...,tj-i,tpf + i,tj + i, tpf}; tj). (95) 
P J P , 

Since k ^ j, the variable set {t\, t k -i, tN+i, tk+i, ■■■,tN} is different from {ti, ...,tj—i,tjf+i,tj+it •••,^jv} only at 
tj and tk- For those Ps such that P{ti, tfe-i, tjv+i, ife+l> ...,ijv} moves tj to the first to positions in the variable 
set {ti, i/v}, the symmetry between t\,...,t m and tjv+i in the function tjv? iiV+i) hi Eq. ((90|) will make 

G(P{t 1 ,...,t k -x,t N+ i,t k +i,...,t N };t k ) = G(P{t 1 ,...,tj^\,t N+ i,tj + i,...,t N };tj). There are totally m(N - 1)1 such 
permutations and they all lead the time integral to 



/ dti...dtN+iG(ti, tj-\)tN+i,tj+\, tj) G*(P'{tx, ij-i, i/v+i, ij+i, — j^jv}; iy)- 

J p, 



(96) 



By Eq.(|48]). it is simply TV. 

For the other permutations that move tj to the position of t m +i, •••,tj\r, it cannot be interchanged with because 
T ^ T s (s = m + l,...,N). Furthermore, by the orthogonal relation in Eq. ([82]) . the time integral is simply zero. 
Therefore, the cross terms are equal to 



/^ 1 ...dt iV+1 ^e i ^-^^(P lfc {t 1 ,..,t iV+1 })g*(P lj -{« 1 ,...,t J v+i}) =m(N-l)Wj2 e 

k^j k^j 



i{S k -Sj) 



(97) 



But because J2 k e = °> 

we have 



E< 



EE 



53 e ' ;5fc E - + !) 

fc 3 

-(N + l). 



(98) 



So the final result is 



P N+1 (lHmV) cx J\f(N + 1)(N - 1)\{N - m) 

= (99) 

For the generalized Hong-Ou-Mandel interferometer, we 
scan the delay of the H-photon relative to the V-photons. 
When it does not overlap with any of the V-photons, 
no interference occurs and Pjv+i is a straight line which 
corresponds to to = in Eq. (f9"9")) with Pjy+i(oo) = (N + 
l)V\f. The value in Eq. (l9"9")) corresponds to the case when 



the delay is zero between the to V-photons and the one 
H-photon and a local maximum interference is achieved. 
So the visibility is 



V N+1 {lHmV) ee 



Pv+i(oo) -P N+1 (lHmV) 



m 
N' 



Pv+i(oo) 



(100) 



Note that this visibility only depends on N and to, i.e., 
the total number N of V-photons and the number to of 
V-photons that overlap with the single H-photon. It is 
independent of the normalization factor J\f or how the 
other N — to photons distribute in time. 

So for a temporal distribution of well separated groups 
of V-photons shown in Fig. 3a, as we scan the location of 
the single H-photon, we will have more dips of various 
visibility (Fig. 3b) and the visibility is m/N when the 
single H-photon overlaps with the group of to V-photons 
that are in one temporal mode and are well separated 
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(a) 



V 



(b) 







FIG. 3: (a) A temporal distribution with well separated 
groups of V-photons and (b) the corresponding normalized 
Pjv+i as the position of the H-photon is scanned. 



and V4(2 x 2) = 1/3. But there are other scenarios like 
2HIV + IV and 1H1V + 1H + IV. We will consider a 
simpler model to include these two scenarios so as to com- 
plete the distinguishability discussion in the four-photon 
case. 

For simplicity, we will again only discuss an un- 
normalized independent four-photon state of the form 



= / duJidijj2duo^dijJ4^{uJi, u>2, u>3, W4) x 

x a\j (uj-l ) a\j ) city (u> 3 ) a\, (w 4 ) 1 0) . ( 10 1) 



from other V-photons. 

In general for a temporal distribution with m partially 
overlapping V-photons, the visibility will be a value less 
than m/N. Therefore, the experimentally measurable 
visibility of the dips can be used to characterize the 
degree of temporal indistinguishability of an N-photon 
state. 



VI. THE GENERAL CASE OF \k„,N v ) WITH 
k > 1 



with ct>2) u>3, W4) in the form of Eq . ([31 



$(wi,W2,W3,w 4 ) = ^i)e MlTl 0(w 2 )e^ T2 x 

x< ? !)( C j3)e' i " 3T3 ( /)(w4)e M4T4 . (102) 



For the NOON state projection measurement with 
V = 4, the field operators at the four detectors are re- 
lated to the input field operators as 



For a more general case of input state of |fcjj, Ny) with 
k > 1, there are many scenarios for the temporal distri- 
bution of the photons. We will start with the four-photon 
case of k = N = 2. 



'Eo(t) = [E H (t)-E v (t)]/2 + .. 
Ej.(t) = [E H (t)+E v (t)}/2 + .. 
E 2 (t) = [E H {t)-iE v {t)]/2 + . 
E 3 (t) = [E H {t)+iE v {t)}/2 + . 



(103) 



A. Four-photon case of \2v,2h) 

This situation was discussed in Ref.[l5j for 4x1 case 
and 2HV x 2HV case. It was shown that V 4 (4 x 1) = 1 



where we omit the vacuum modes. The four-photon de- 
tection probability at the four detectors is related to the 
following correlation function: 



r< 4 >(fi,f 2 ,*3,*4) = (^(t4)^i t (t3)S 2 t (i2)S 3 t (il)^3(il)^2(t2)^l(t3)S (< 4 )). (104) 

Again, it is easy to first calculate -^(ii)^^)^^)-^^)!^)- For this, we expand E 3 (ti)E2(t2)Ei(t3)E (t4 : ): 

E3(ti)E2(t2)Ei(t 3 )E (t4.) = [{VVHH - HHVV) + i(VHVH + HVHV) - i(HVVH + VHHV)]/W, (105) 



where H = E H ,V = E v 
HHVV\<S>i) 


and we keep the time ordering. For the state |<E> 4 ) in 

= [G(tl,t2,t3,t4,) + G(t2,tl,t3,ti) + G(ti,t 2 ,t4,t 3 ) + 


Eq.(jl()l|), we have 
GfaMMMW), 


(106) 


VVHH\$ A ) 


= [G(*3,t 4 ,tl,t 2 H 


- G(t 4 , t 3 , ti,t 2 ) - 


r-[G(t 3 ,t4,*2,*l)4 


G(t 4 ,t 3 ,t 2 ,t 1 )}\0), 


(107) 


HVHV\1> 4 ) 


= [G(tx,t3, t 2 ,t 4 ) - 


- G(t 3 , tx,t2, t 4 ) - 


\-G{t 1 ,t 3 ,t A ,t 2 ) + 


G(h,t 1 ,t 4 ,t 2 )}\0), 


(108) 


VHVH\<S> 4 ) 


= [G(t2, ti, ti, t 3 ) - 


- G(t2,ti, t 3 , t\) - 


{-G(t A ,t 2 ,h,t3) + 


G(t 4l t 2 ,t 3 ,t 1 )}\0), 


(109) 
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HVVH\$ 4 ) = [G(t 1 ,t 4 ,t2,t3) + G(h,U,t3,t2) + G(t A ,h,t 2 ,t 3 ) + G{t 4 ,tx,t 3 ,t 2 



(110) 



VHHV\<5>4 



[G(t2,t 3 ,t 1 ,U) + G(t2,t3,U,t 1 ) + G(t 3 ,t2,h,U) + G(t 3 ,t 2 ,U,t 1 )}\0) 



with 



G(t X ,t2,t 3 ,t 4 ) 



I 



(111) 



(112) 



For the <I>-function given in Eq. (|102p . the above G- 
function is simply 

G(ti, t2, t 3 , t 4 ) 

= g(h ~ TMh ~ T 2 )g(t 3 - T 3 )g(t 4 - T 4 ). (113) 

Four-photon coincidence probability is proportional to a 
time integral of the correlation function r' 4 ^: 

Pa = J dt^dhdUT^^MMM)- (H4) 

Next, we will evaluate P4 for various scenarios of 
I 



photon distinguishability. To describe the four scenar- 
ios discussed in the beginning of this section, we in- 
troduce three delay parameters: AT, ATy, ATh so that 
T 2 = T 1 + AT H , T 3 = T 1 + AT, T 4 = T 3 + ATy. There- 
fore, AT is for the delay between the H-photons and the 
V-photons and AT H (v) for the delay between the two 
H(V)-photons. When AT = ±00, there is no overlap be- 
tween the H- and V-photons and no interference occurs. 
This sets up the baseline for evaluating the visibility of 
interference. We start with the 4x1 case: 

(i) ATh = = ATy. There is an exchange symmetry 
between ti,t2 and between t 3 ,t 4 in G(t x ,t 2 ,t 3 ,t 4 ) with 



G(t x ,t 2 ,t 3 ,t 4 ) = g(t x - T x )g(t 2 - T x )g(t 3 -T x - AT)g(t 4 -T x - AT). 



(115) 



So we have 
E 3 



{h)E 2 {t2)Ei (t 3 )E Q (t 4 )\<f> 4 ) = j{[G(t x , t 2 , t 3 , U) - G(t 3 ,t 4 , t x ,t 2 )} + i[G(t x ,t 3 ,t 2 ,h) + G(t 2 ,U, h,t 3 )] 

-i[G(ti,U, h, h) + G(t 3 , ta, ti, **)]}|0). (116) 



After the time integral, we obtain 



4 2) (r)= ( / dw|^( W )| 2 e 



1 



(120) 



P 4 (AT) = - [3^4 — ££4, (AT) + £ 4 (AT)] (117) 

o 



with 



Note that £i 1} (0) = £ y 4 '{0) = A 4 and £ y 4 '(oo) = 

(2) 

Ef'(oc) = 0. As we scan the relative delay AT between 
(118) the H- and V-photons, the four-photon coincidence will 
show an interference dip all the way to zero when AT = 
0, which corresponds to the case of T x — T 2 = T 3 = T 4 
or the 4x1 case. So the visibility is 100% for the 4x1 



duj\(j>(uj)\ : 



f| 1) (r)=( I dcj\^e 



duo\4>{uj)\ z , (119) case. 



(ii) ATh = but ATy >> 1/Aw. In this case, the two V-photons are well separated and we have 
G{h,t 2 ,t 3 ,U) = g(t x - T x )g{t 2 ~ T x )g(t 3 -T x - AT)g(t 4 -T x - AT V - AT). 



(121) 



When AT = 0, there is an exchange symmetry between {t x ,t2,t 3 } in G(t x ,t 2 ,t 3 ,t 4 ). This is the 2H1V + IV case. 
But for arbitrary AT, there is only a permutation symmetry between t x ,t 2 in G(t x , t^, t 3 , £4). Then we have 



E 3 (t x )E2(t2)E x (t 3 )E (U)\$4 
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g t 2 , t 3 , t 4 ) + G(tl, t2,t 4 , £3) — G(^3, £ 4 , ti, <2) — G(t 3 , ti, *2, ^l)] 

+i[G(ti,ta, * a , **) + *3, *4, t 2 ) + G(i 2 , < 4 , h,t 3 ) + G(t 2 ,t 4 , t 3 ,ti)] 

-i[G{t u t 4 , t 3 , h) + G(t 3 ,t 2 , t lt t A , ) + G{t lt t A , t 2 ,t 3 ) + G(t 3 ,t 2 , u, h, )]}|0). 

I 



(122) 



When AT = ±00, there is no overlap between all the 
terms in Eq. (|122p so that all the cross terms are zero 
after the time integral in Eq. (|114p . So we have 



P 4 (AT = ±00) = 3A/16. 



(123) 



On the other hand, when AT = 0, there is an ex- 
change symmetry between {£1,^2,^3} in G(ti,t 2 ,t 3 ,t 4 ). 
So Eq. (fl%|) becomes 

£ 3 (ii)£2(t2)£i(t3)£o(i4)|$4) 

G(ti,t 2 , t 3 , t 4 ) + G(ti,t 2 ,t 4 , t 3 ) 

-G(t 3 , U, h,t 2 ) - G(t 3 , U, h, ti)] |0) (124) 

and there is no overlap between all four terms above. 
After the time integral, we obtain 



P 4 (AT = 0) = A/16. 
So the visibility is 

V 4 {2H1V + IV) = 2/3. 



(125) 



(126) 



for the 2H1V+1V case. In fact, there is another 2H1V+ 
IV case when the two H-photons overlaps with the other 
V-photon and AT = — ATy. In this case, we have the 
exchange symmetry between {t\, t 2 , t 4 } in G{t\,t 2 , t 3 , t 4 ) 
so that 



P 4 (AT = -ATy) = At/16, 



(127) 



which also gives V 4 (2P1V + IV) = 2/3. 

{in) AT H = ATy = T » 1/Au. This is the 
1H1V + 1H1V or the 2 x 2 case when AT = and 
we have the exchange symmetry between {t\, t 3 } and be- 
tween {t 2 ,t 4 }. But for AT = ±00, there is no overlap 
between any two of the 24 terms in Eqs. (|106Hlll|) . So 
we have after the time integral in Eq. (|114[) 

P 4 (AT = ±00) = 24A/16 2 = 3A/32. (128) 

When AT = 0, on the other hand, we have 



E 3 {ti)E 2 {t 2 )E l {t 3 )E {U)\^ i ) 



G(ti,ta, U, ta) + G{t 3 ,h,t 2 , t 4 ) - G{t u t 4 , t 3 , t 2 ) - G(t 4 , h,t 2 ,t 3 ) |0). (129) 
I 



The above four terms have no overlap so that we obtain 

P 4 (AT = 0) = 4A/8 2 = A/16. (130) 

Therefore, the visibility for the 2x2 case is simply 

V 4 (2 x 2) = 1/3. (131) 

(iv) \ATh — ATy I >> 1/Aw and |AT H |,|ATy| >> 
I/Aw. As we scan AT, there is an exchange symme- 
try only in one pair of the variables between {£1,^2} 
and {t 3 ,t 4 }, that is, between {ti,t 3 } when AT = 0, 
or between {ii,i 4 } when AT = —ATy, or between 
{t 2 ,t 3 } when AT = ATh, or between {t 2 ,t 4 } when 
AT = AT H - ATy. This is the {I HIV + IH + IV) 
case. In all these cases, 8 out of 24 terms in Eqs. (|106t - 
II 1 ID are cancelled in Eq. (|105p and the remaining ones 
are orthogonal to each other so that we have 



The situation when AT = ±00 is same as Eq. (|128l) . 
Therefore the visibility is 



V 4 (1P1V + 1H+ IV) = 1/3. 



(133) 



P 4 (AT = 0) = 16A/16 2 = A/16. 



(132) 



for the 1H1V + 1H + IV case. 

These are all likely distinct scenarios. We summarize 
the visibility in Table I. Although visibility is derived 
with a specific ^-function in Eq. (|102p . in general, visibil- 
ity is the same regardless of the form of <& as long as it is 
such that G(t\, t 2 ,t 3 , t 4 ) has the required exchange sym- 
metry in each scenario listed above. The intermediate 
situations will not have any symmetry in G(t±, t 2l t 3l t 4 ) 
and thus have very complicated dependence on the vari- 
ous permutations of u>2, k>3, u> 4 ). Ref.[15| discussed 
the intermediate scenario from the 2x2 case to the 4x1 
case. Indeed, the visibility depends on the quantity £ /A, 
which defines the degree of pair distinguishability. Xiang 
et al [12] realized the 2x2 and the 4x1 cases experi- 
mentally and confirmed the visibility in Table I. 
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TABLE I: Visibility for 2 H-photons and 2 V-photons input 

2H2V 2H1V+1V 1HV+1HV 1HV+H+V 

1 2/3 1/3 1/3 



TABLE II: Visibility for 2 H-photons and 3 V-photons input 

2H3V 2H2V 2H1V 1H3V 1H2V 1H2V HV+V HV+V 
+V +2V +H +HV +H+V +HV +H+V 

1 5/6 1/2 3/4 5/12 1/2 1/3 1/4 



B. The Special Cases of \2 H ,S V ), \2 H Av), and 
|3ir, 3v) 

Following the same line of derivation but in a much 
more complicated fashion, we may find the visibility for 
all the scenarios for the input states of |2# , 3y), \2h, 4y), 
and \3h , 3y). We list the likely scenarios below and tab- 
ulate the visibility for each scenarios in Tables II-IV. 

1. The Case of |2 H ,3v) 

The case of |2#,3y) has 8 different scenarios. They 
are 

(i) \2H3V),\2H2V +IV),\2HIV + 2V), and 

(ii) \IH3V + IH), \IH2V +IH1V), \IH2V + 1H + IV), 
\1H1V + 1H1V + IV), \1H1V + 1H + 2V). 

Their visibilities arc listed in Table II. 



2. The Case of\2 H Av) 

The case of |2#-,4y) has 12 different scenarios. They 
are 

(i) \2HW), \2H3V + IV), \2H2V + 2V), \2H1V + 3V), 
and 

(ii) \1H4V + 1H),\1H3V + 1HIV), \IH3V + IH + 
IV), \1H2V + 1H2V), \1H2V + IHIV + IV), \IH2V + 
IH + 2V), \IHIV + IHIV + 2V), \1HIV + IH + 3V). 

The scenarios with different visibility are listed in Ta- 
ble III. \1H2V + IHIV + IV) and \IH2V + IH + 2V) 
have the same visibility of 2/5 as \1H2V + 1H2V). 

In general, they follow the trend that smaller visibil- 
ity corresponds to less photon overlapping. However, 



there are exceptions: 1H2V + HV has less visibility than 
1H2V + IH + V in Table II and 1H3V + HV has less 
visibility than IH3V+1H+1V in Table III. So the run- 
away HV does not help when H and V overlap in these 
cases. 

3. The Case of\3 H ,3 V ) 

There are totally 11 different scenarios in the special 
case of \3h, 3v)'- 

(i) \3H3V),\3H2V + V),\3H1V + 2V); 

(ii) \2H2V+1HIV), \2H2V+1H+IV), \2H1V+1H2V), 
\2H1V + IHIV + IV), \2H1V + 1H+ 2V); 

(iii) \IHIV +IHIV +IHIV) , \IHIV +IHIV +IH +IV) , 
\IHIV + 1H+ IV + IH + IV). 

In Table IV, we list the visibility for most of the sce- 
narios. \2HIV + IHIV + IV) and \2HIV + IH + 2V) 
have the same visibility of 2/5 as \2H1V + 1H2V) and 
are not listed. As can be seen, anomaly occurs for 
\2H2V+1H1V) and \2H2V+1H+1V) where visibility 
is bigger for the case with less photon overlap. The sce- 
narios of \3H3V), \2H2V + IHIV), and |3 x HV) were 
observed experimentally by Xiang et al. [22j with the 
corresponding visibility in Table IV. 



C. General Formula for the Visibility 

The most general case is when the input state is in 
the form of \kn,Ny) with k < N. The most general 
scenario is when the k H-photons don't overlap in time 
but rather are split into r temporally well separated sub- 
groups with kj indistinguishable photons in the jth group 
and kx + ... + k r = k. We also divide the N V-photons 
into r + 1 subgroups with rrij V-photons overlap in time 
with the jth. H-photon group. The rest N — mi — ... — m r 
V-photons are in a separate group by themselves. The 
wave function for these N + k photons will satisfy the 
permutation symmetry relation similar to Eq. (|26[) for the 
overlapping photons and the orthogonal relation similar 
to Eq. (|2~T!) for the well separated photons. 

The derivation of the general formula for the visibility 
in the (N + fc)-photon NOON-state projection measure- 
ment is very complicated and lengthy. It follows the gen- 
eral line of argument as that leading to Eq. (|100[) . We will 
present the detailed procedure elsewhere [23| but only 
give the result as 



Ec- 1 )'" 1 
i=i 



i 

E 



p. 



/~ti T (il) (i r ) 

.C k r m\ ...m r 



h\..A r \ ) (N + k- 1)...(N + k~iy 



(134) 
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TABLE III: Visibility for 2 H-photons and 4 V-photons input 



2H4V 


2H3V 


2H2V 


2H1V 


1H4V 


1H3V 


1H3V 


1H2V 


2xHV 


1H1V 




+V 


+2V 


+3V 


+H 


+HV 


+H+V 


+1H2V 


+2V 


+1H+3V 



1 9/10 7/10 2/5 4/5 1/2 3/5 2/5 3/10 1/5 



TABLE IV: Visibility for 3 H-photons and 3 V-photons input 



3H3V 


3H2V 


3H1V 


2H2V 


2H2V 


2H1V 


HVx3 


HVx2 


HV+V 




+V 


+2V 


+HV 


+H+V 


+1H2V 




+H+V 


+H+H+V 


1 


9/10 


3/5 


3/5 


7/10 


2/5 


2/5 


3/10 


1/5 



where = = m^ m \ m^' = m(m — l)...(m — i + 

1), and C^ 1 = (N + M)\/N\M\. For the special case 
of k = 1, Eg. (|134p recovers the expression in Eq. (|100l) . 
Furthermore, we can easily check that the formula in 
Eq. (|134p indeed leads to the visibility values in Tables 
I-IV. 



VII. CONCLUSION AND DISCUSSION 

The complementary principle of quantum interference 
is demonstrated in a quantitative way in multi-photon in- 
terference where photons can be categorized by their tem- 
poral distinguishability. The temporal indistinguishabil- 
ity of photons in turn can be characterized by the permu- 
tation symmetry in the multi-photon wave function while 



the temporal distinguishability by the orthogonality of 
the permuted wave functions. Generalization to other de- 
grees of freedom such as spatial modes is straightforward. 
Although the above conclusions were made on photons, 
they should apply to any bosons as well as fermions so 
long as the occupation number of each mode is less than 
or equal to one. 
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